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Abstract 

Dye sensitised nanocrystalline solar cells (Gratzel cells) 
have achieved solar-to-electrical energy conversion efficien- 
cies of 12% in diffuse daylight. The cell is based on a 
thin film of dye-sensitised nanocrystalline Ti02 interpen- 
etrated by a redox electrolyte. The high surface area of 
the TiC>2 and the spectral characteristics of the dye allow 
the device to harvest 46% of the solar energy flux. One 
of the puzzling features of dye-sensitised nano-crystallinc 
solar cells is the slow electron transport in the titanium 
dioxide phase. The available experimental evidence as well 
as theoretical considerations suggest that the driving force 
for electron collection at the substrate contact arises pri- 
marily from the concentration gradient, ie the contribu- 
tion of drift is negligible. The transport of electrons has 
been characterised by small amplitude pulse or intensity 
modulated illumination. Here, we show how the transport 
of electrons in the Gratzel cell can be described quan- 
titatively using trap distributions obtained from a novel 
charge extraction method with a one-dimensional model 
based on solving the continuity equation for the electron 
density For the first time in such a model, a back reaction 
with the ±3 ions in the electrolyte that is second order in 
the electron density has been included. 

PACS: 72.20.Jv,72.40+w,73.50.Pz,73.63Rt 
Keywords: photovoltaics, nanorystalline, transport. 



1 Introduction 

Dye-sensitised nanocrystalline cells were developed by 
Gratzel and coworkers jy, gj. Their operation is shown 
in figure 0. Solar radiation is absorbed by a monolayer 
of a ruthenium based dye adsorbed on the surface of a 
porous nanocrystalline film consisting of electrically con- 
nected particles of ~ nm in diameter of Ti02- On excita- 
tion of the dye by the radiation, electrons are transferred 
from the HOMO to the LUMO band of the dye, the lat- 
ter lying just above the conduction band of the Ti02, and 
then injected into the Ti02 on timescales thought to be 
~ fs B El. The Ti02 film is permeated by an electrolyte 
which contains iodide ions (I - ) and iodine in the form 
Ig~ . The dye is regenerated by electron transfer from the 
I - ions, leading to the formation of 1 3". In turn, the Ig~ 
ions are regenerated by electron transfer from the cath- 
ode, facilitated by a layer of platinum which acts as a 
catalyst. Electrons injected into the Ti02 diffuse through 
the nanocrystalline film to the anode. 

Whilst diffusion is slow in the electrolyte, the high con- 
centrations of Ig" and I - ions result in loss free hole trans- 
port. Hence much of the research on these cells has fo- 
cussed on electron transport in the nanocrystalline film 
since injected electrons are slowed down by trapping at the 
surface of the Ti02 particles and may back react through 



recombination with Ig" ions. As the surface area/volume 
ratio is so high, it seems probable that traps at or near the 
surface dominate over those in the bulk. Back reaction 
with the dye ions occurs both through geminate reactions 
and from traps. It is likely that back reactions with the 
dye ions are reduced by efficient scavenging of the photo- 
generated holes by the I - ions 0, so back reaction with 
I3 ions will be a significant cause of loss of efficiency. Mea- 
surements of the decay of the total electron concentration 
at open circuit show that the rate of back reaction with 
17 ions is second order in the total electron concentration 
p, pfl , and possible back reaction sequences consistent with 
this observation were proposed in M . 

The response of the cell to time varying light intensities 
can yield much information about the density of traps, 
their energy distribution and the rate and nature of the 
back reaction. This can be done through the large ampli- 
tude measurements of electron decay in (3, |§], but these 
have the disadvantage that the whole trap distribution is 
sampled in a single measurement. More information can 
be obtained by employing intensity modulated photocur- 
rent spectroscopy (IMPS) and intensity modulated photo- 
voltage spectroscopy (IMVS) which involve superimposing 
a small modulated component of light intensity on a dc 
component I |2| || . A plot of the imaginary part of the 
modulated photocurrent against the real part of the same 
quantity gives a semicircle in the negative imaginary half 
of the Argand diagram. The imaginary component is neg- 
ative as the response lags behind the light intensity and is 
the majority carrier pj. The magnitude of the imaginary 
response is a maximum when the cell is being perturbed at 
a frequency u> m i n which is the inverse of its response time 
t. Here, r is the mean transit time for electrons through 
the Ti02 network, allowing for electron scattering which 
is included in the bare diffusion coefficient Db are (ie with- 
out traps), and multiple trapping. A similar result can be 
derived for the photovoltage. 

We have used this approach to find how r from IMPS 
spectra varies with the background intensity Jo- Increasing 
Iq results in an increase in the conduction band density 
and hence the electron quasi- Fermi level Ep n . As Ep n 
resides in the bandgap, changing 1$ scans Ep n through the 
trap levels. Given that the response is dominated by states 
near Ep n , IMPS and IMVS provide a means of probing 
the trap density of states. This was shown || for a first 
order backreaction with the assumption that the gradient 
of the electron quasi- Fermi level is constant. 

Below, we show how uj m i n can be calculated from a 
multiple trapping model in which electron transport is 
mediated by the conduction band and is interrupted by 
trapping. An alternative approach is the hopping model 
in which it is assumed that the electrons are always in a 
trapped state and transport is achieved by hopping from 
one trap to another. Recently Nelson et al |10| looked at 
both models using a random walk Monte Carlo model util- 
ising the continuous time random walk method of Scher 
and Montroll 1 1 1 where steps to all nearest neighbours are 
equally likely and the waiting time depends only on the 
energy of the initial site, and which included the back re- 



action with I^~ ions. They showed that the multiple trap- 
ping model is more capable than the hopping model of 
explaining experimental data on charge decay due to re- 
combination. 

To calculate the IMPS spectra, we solved the continuity 
equations for the conduction electron density and trapping 
probability. Our approach is a generalisation of that of 
Vanmaekelbergh and de Jongh || where we have included 
back reactions of first and second order and avoided the 
assumption that the gradient of the electron quasi-Fermi 
level is constant. 

A common treatment of the transient response, used 
is to include the traps through an effective dif- 
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fusion coefficient which depends on the conduction elec- 
tron density or equivalently light intensity. For the diffu- 
sion coefficient Db are we have taken the value measured 
on nanocrystalline anatase TiC>2 which has a small defect 
concentration. This is because we have explicitly allowed 
for the influence of electrolyte and defects in our model so 
they should not be also included in the diffusion coefficient. 
In , expressions for an effective diffusion coefficient and 
effective transition time that explicitly show the influence 
of the traps and back reaction were obtained from a sim- 
ilar model to ours but for a single trap level and a first 
order back reaction. However, in our model the physical 
meaning of these quantities is not clear. 

2 Method 

The rate of change of electron density n(x, t), where x is 
the distance from the anode and t the time, is determined 
from the continuity equation for conduction electrons: 
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Here the first term on the rhs is the photocurrent due to 
with diffusion coefficient without traps D oare , the second 
term is the generation rate from light of intensity 1$ with 
an absorption coefficient a, the third term is the back 
reaction rate for conduction electrons with the IJ ions at 
a rate of fc c f, M , where fj, = 1 for 1st order recombination, \x 
= 2 for 2nd order recombination and ndark the conduction 
electron density in the dark, and the fourth term is the net 
trapping rate whose terms are defined below. Equation pi is 
coupled with the continuity equation for trapped electrons. 

df 
Nto < -q7 >=< hra P n(l - f) - k detrap N t of > 

— ktbn{n trap — n trapdark) W 

where Nto is the density of traps, / the probabilty that 
a trap level is occupied, k tra p the trapping rate, kdetrap 
the detrapping rate and fc t b M the back reaction rate for 
trapped electrons with the Ij ions. The angled brackets 
<> represent an average over the trap distribution, for 
any quantity A: 
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where s(Et)&Et is the probability of a trap site existing 
in the energy range [Ef, Et+AEt] in the bandgap E g , 
and is normalised so that it integrates to unity over the 
bandgap. Thus, the trapped electron density 
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and ntrapdark is n tra p in the dark. Detrapping by thermal 
excitation is assumed, so 
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where ks is the Boltzmann constant. 

An exponential trap distribution was used here in com- 
mon with other authors ( p0[ ^3|), using the parameter (3 
to represent the width of the distribution in units of the 
thermal energy ksT : 
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If (3 is large, say 0.3, most of the traps sit up near the 
conduction band, so detrapping will take place on average 
relatively quickly. As (3 is reduced, the probability of hav- 
ing deeper trap levels increases corresponding to longer de- 
trapping times. The limits — > oo and [3 — ► correspond 
respectively to all traps situated just below the conduction 
band and to a uniform trap distribution throughout the 
bandgap. We have also looked at Gaussian distributions 
as recent transient current measurements E3] suggest that 
the trap distribution may be peaked but these results will 



be published separately |15 
The photocurrent density 
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where q is the magnitude of the electron charge. At the 
anode, 

j P hoto(0) = qk ext {n(0) - n dark ) , (8) 



where k ext = ^/fcsT/(2TO*7r) for IMPS spectra (short cir- 
cuit) and k ex t — for IMVS spectra (open circuit). This 
is the Schottky boundary condition |p_6|| with a recombi- 
nation velocity v r = k ext . At the cathode side of the TiC>2 
network, 

jphatoid) = . (9) 

The incident photon to electron conversion efficiency 
IPCE is j p hotoss(x = 0)/qI a where jphotoss is the steady 
state photocurrent density. If trapping effects are ne- 
glected, it can be shown from equation [j] that the IPCE 
is equal to 1 — exp(— ad). The ac photocurrent conversion 
efficiency is $(w) = j p hoto{x = Q)/qSIo, We solved equa- 
tions [y and g in the steady state and for the modulated 
intensity SIq exp(iwt) (where SIq ~ O.Ol/o so that only the 
term linear in SI was included). Hence $(w) is complex. 

As the electrostatic potential is screened out, classical 
statistics gives 



n = N c exjp(E Fn /k B T) 



(10) 



Table 1: Parameter values 



Symbol 


Value 


-Ubare 


1.0 x 10" b nr^s" 1 


< 


9m 
6.775x 10 26 itT 3 


a 


2.3 x 10 5 m- 1 


T 


300 K 


J^Fndark 

fctrap 

d 


-1.0 eV 

1.0 x 10 7 s" 1 
6 x 10~ 6 m 


k c b2 = k t b2 


2 x 10- 25 m 3 s _1 
0.2 



Table 2: Values n(x = 0) at different background intensi- 
ties Jq (m~ 2 s _1 ) and trap densities N t o(m~ 3 ) 



where the conduction band density of states 



N r . 



m* n k B T 



2-kTi 



3/2 



(11) 



where m* is the effective mass and Ti is Planck's constant 
divided by 2it. 

To solve the equations, a relaxation method for a two 
point boundary value problem [[17] was adopted. Sec- 
ond order kinetics resulted in nonlinear functional equa- 
tions for the steady state and modulated components of 
/, which were solved using the Newton-Raphson method 

The parameter values are given in Table Ql. Di> are and 
to* (in units of the bare electron mass too) were taken 
from ref [ [18[ . We deduced a by assuming an IPCE of 
0.75 IT9). Our dark Fermi level EFndark comes from an 
estimate of 10 11 m~ 3 for ndark M- For an estimate for 
ktrap, we assumed that the trap cross-section a would be 
in the region of 2xl0 _22 m -2 , and that k tra p = Ntoavth 
where the thermal velocity of the electrons Vth = V^TTVr, 
but to keep the model simple, we made ktrap independent 
of N t Q. Our chosen values for k c b2,ktb2 and j3 are discussed 
at the beginning of section |3[ 

3 Results 

Our results are shown for a limited range of values of the 
trap density N t0 , trap distribution width and back re- 
action rate k t b2- We have set k t b2=k c b2 on the grounds of 
simplicity. Our results are insensitive to k c j,2 as most of 
the electrons are trapped. We have chosen a value for ktbi 
100 times less than the value deduced from experiment H 
as this value was found to give better agreement between 
our IMPS predictions and experimental results. 

For the parameter values shown, our calculations of 
^min give the best agreement with experiment for a range 
of values. How our predicted spectra vary with j3 and 
ktb2 will be presented in a separate paper [g_5|, where our 
results for the IPCE as a function of intensity will also 
be given. We have not attempted to find a set of values 
of N t Q,ktb2 an( i /? to give the closest fit to experimental 
data as the model does not seem accurate enough to jus- 
tify this, and it is more instructive to see the trends in 
ijJmin as Nto,ktb2 and /? are varied. We chose the values of 



Nta 


7o=lxl0 15 


7o=lxl0 19 


2.3xl0 23 


9.41 xlO 10 


8.35xl0 14 


2.3xl0 24 


8.53xl0 10 


8.34xl0 14 



N t Q from noting that the density of states obtained from 
charge extraction measurements is consistent with an 
exponential distribution of the type given in equation o 
for which N t0 = 2.3xl0 24 uT 3 and /3=0.14 |(j). However, 
more recent data suggests that j3 may be as high as 0.25 
pl[ . Our chosen value for (3 lies within this range, and it is 
at the bottom of the range of values suggested by Nelson 
et al liri] when modelling the decay of photoinduced dye 
cations. 

3.1 Steady State Results 

These results were obtained by setting the lhs of equations 
pi and to zero. In figure @, we show an example profile 
for the conduction band density n(x). This figure demon- 
strates the large density gradient that drives the electrons 
to the anode. The photocurrent is determined by the trap 
distribution at the anode where the electron density is 
much lower compared with its value on the cathode side. 
This means that there will be significant differences com- 
pared to IMVS spectra at the same intensity, for which 
n(x) is almost constant throughout the Ti02 film because 
of the boundary conditions forcing the gradient of n(x) 
to be zero on both sides of the film. As the profiles for 
n(x) are the same shape at all intensities and n(x = 0) 
is the only value needed for the photocurrent and IMPS 
results, we have only given our results for 77,(2; — 0) which 
are shown in Table g. These results show that n(x = 0) is 
nearly proportional to the light intensity, so that the IPCE 
is almost constant as would be expected for the low back 
reaction chosen. It is interesting to note that n(x = 0) is 
insensitive to the trap density at both light intensities. 

The difference in behaviour for the different trap densi- 
ties is rather seen in the steady state values for / at the 
anode which are shown in figure [| for the same trap den- 
sities and light intensities as in Table |[ This graph shows 
the expected Fermi-Dirac distribution at the lower inten- 
sity, although for the lower trap density the spread in / 
is greater than for the higher trap density. This leads to 
a larger nt ra p and hence back reaction rate compared to 
the solution for the higher trap density. From equation 0, 
setting the left hand side to zero in the steady state, it can 
be seen that a large back reaction forces / to be less than 
one at low trap energies, and this effect increases, the more 
ntrap exceeds n tra pdark- The same phenomenon is seen at 
the higher trap density but here the difference between 
ntrap and n tra pdark is much greater at this density, forcing 
the observed low value of / at low trap energies. 

The plateau between trap energies of -0.7 and -0.9 eV at 
the higher intensity and trap density resembles closely the 
plateaux predicted by Simmons and Taylor [2^] in calcu- 
lations of / in the steady state for bipolar semiconductors 



with traps. Here, we have the same system except that 
it is unipolar and there is a second order back reaction, 
so the fact that we see the same features is reassuring. 
What is different from reference |2^] is our result that / is 
less than unity even at trap levels deep into the conduc- 
tion band. This is a direct consequence of the assumption 
that the back reaction rate is independent of the trap en- 
ergy Et, and may therefore be an artifact of the model. 
However, it does not unduly affect the results as what is 
important is the overlap with the trap distribution s(Et) 
and this decreases rapidly for low values of Et- 



3.2 IMPS Results 

Figure [| compares the predicted IMPS spectra for two 
light intensities. The spectrum is approximately semicir- 
cular. Im($(w ra i„)) at the higher intensity is about three 
times its value at the lower intensity. This can be qualita- 
tively explained as being due to the greater losses due to 
the back reaction at the higher intensity. 

The semicircle is a result of a single value for r even 
though there is a distribution of trap levels; Vanmaekel- 
bergh et al || note that this is a result of the capture rate 
for a given trap level varying inversely with the residence 
time. However at the higher intensity, where the back re- 
action is more important, the semicircle is more distorted 
as seen from figure |[ This distortion suggests that there 
is a dispersion in the electron arrival times at the anode 
H . For each spectrum, it can be seen that there is a single 
value of lOmin within the range of frequencies calculated. 
Whilst another value exists at much higher frequencies, 
this has not been calculated as it would occur well outside 
the frequencies accessible to experiment. 

Our main result is shown in figure ra where our predic- 
tions for Lo m in are compared with experimental data from 
p3| . The magnitude of ui m in = 1/r is determined by the 
time taken in trapping and detrapping processes. As I 
increases, the traps fill up (see figure 0) so that detrap- 
ping occurs from the shallower trap levels at which it takes 
place faster from equation |l From figure 0, it is clear that 
at any given intensity over the range considered, u> m i n de- 
creases as iVto increases. This is because if there are more 
traps, the likelihood of an electron being trapped and then 
detrapped en route to the anode increases, slowing it down 
so increasing r. 

The good agreement with experiment is encouraging. 
In our calculations for all three trap densities, we have 
achieved the same gradient of approx 0.8 for log 10 (tL> m i„) 
vs log 10 (Jo) as in the experiments at the higher intensities. 
However, unlike our calculations, the experimental results 
flatten off at the lower intensities. This effect seen in our 
calculations for larger values of k t b2 , although in our calcu- 
lations the flattening off exceeds that seen experimentally. 
The flattening occurs because the trap probability in the 
steady state / does not differ greatly from its dark value 
except at the higher intensities. 



4 Conclusions 

We have shown for the first time with a one-dimensional 
model including explicitly trapping and back reaction with 
Ig ions that the assumption of a back reaction which is 
second order in the electron density gives values of to m in 
deduced from IMPS spectra which are in good agreement 
with experiment. Our model will be useful in obtaining 
quantitative information from intensity modulated spec- 
tra on the trap distribution and back reaction and hence 
establish how different methods of cell fabrication yield 
different conversion efficiencies. Many questions however 
remain unanswered, for example whether the same param- 
eters will give equally good agreement with IMVS spectra 
and whether a suitable set of parameters could give good 
agreement with experiment if a first order back reaction 
were assumed. 
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Figure 2: Conduction band density profile n{x) where 
x is the distance from the anode for an intensity Iq of 
lxl0 19 m- 2 s _1 and a trap density N t0 of 2.3xl0 24 uT 3 
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Figure 1: Schematic diagram of solar cell 
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Figure 3: Trap occupation probability / as a function of 
the trap energy relative to the conduction band energy Et 
for intensities Iq of 1 x 10 15 m _2 s _1 and 1 x 10 19 m _2 s _1 and 
trap densities N t0 of 2.3xl0 23 m" 3 and 2.3xl0 24 m" 3 
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Figure 4: Imaginary part vs real part of the ac pho- 
tocurrcnt conversion efficiency $ at intensities Iq of 
lxl0 15 m _2 s _1 (squares) and lxI0 19 m _2 s _1 and a trap 
density N t0 of 2.3xf0 24 m" 3 (circles) 



10 d 



10" 



: experiment 
'N t0 =2.3x10 23 m" 3 




10' 



10" 



10" 
Intensity l m~V 1 



10" 



10 



Figure 5: Minimum frequency to m in vs background inten- 
sity Iq for trap densities N t o of 2.3xl0 23 m -3 (diamonds), 
7.3xl0 23 m~ 3 (squares) and 2.3xl0 24 m~ 3 (circles) com- 
pared with experiment (crosses) 



